Disorder-induced freezing of dynamical spin fluctuations in underdoped cuprates 
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We study the dynamical spin susceptibility of a correlated d-wave superconductor (dSC) in the 
presence of disorder, using an unrestricted Hartree-Fock approach. This model provides a concrete 
realization of the notion that disorder slows down spin fluctuations, which eventually "freeze out" . 
The evolution of disorder-induced spectral weight transfer agrees qualitatively with experimental 
observations on underdoped cuprate superconductors. For sufficiently large disorder concentrations, 
static spin density wave (SDW) order is created when droplets of magnetism nucleated by impurities 
overlap. We also study the disordered stripe state coexisting with a dSC and compare its magnetic 
fluctuation spectrum to that of the disorder-generated SDW phase. 

PACS numbers: 74.25.Ha,74.72.-h,75.10.Jm,75.40.Gb 



Introduction. Recent studies of layered high-Tc su- 
perconductors have provided evidence for nematic order 
of the underlying electronic state[l|-|5|- This develop- 
ment has reinvigorated studies of stripe phases in the 
cuprate materials, and shown the important role of dis- 
order. Here we focus on the magnetic excitations within a 
model that allows explicitly for disorder-inducedpinning 
of incommensurate (IC) stripe spin fluctuations [6| • A de- 
tailed understanding of the magnetic fluctuations in high- 
Tc materials is likely to be important for understanding 
the mechanism of superconductivity itself. Remarkable 
experimental progress in recent years has provided a con- 
verging picture of an hourglass- like dispersion of IC mag- 
netic fluctuations which appear to be universal except for 
a low-energy branch which depends, e.g, on doping and 
material-specific parameters [TJ. All materials appear to 
exhibit a spin-gap beyond some critical hole concentra- 
tion Xc- By contrast, in the underdoped regime, the fluc- 
tuations soften and eventually freeze out to form a glassy 
spin state. This doping dependence seems to apply to 
both "clean" cuprates like YBa2Cu306-j-x (YBCO) where 
quasi-static SDW order is foundlsj, IsMllI]. and intrinsi- 
cally disordered materials like La2-a;Sr2:Cu02 (LSCO) 
where the static spin correlations are long-range, and per- 
sist for a large doping range well into the dSC dome|12|. 
The spin gap in the clean materials appears to persist to 
substantially lower doping, however. 

Experimentally, it is known that an applied magnetic 
field or impurity substitution cause additional slowing- 
down of spin fluctuations |14| - |2C)| . For example, in Ref. 
18l it was shown how substitution of nonmagnetic Zn ions 
for Cu in 15% doped LSCO shifted spectral weight into 
the spin gap, and eventually, for enough Zn (^ 2%), gen- 
erated elastic magnetic peaks in the neutron response. 
In YBa2(Cui_j,Zny)306.97 with y = 2% evidence for a 
similar in-gap Zn-induced mode was observed by Sidis 
et aZ.[13| Upon increased temperature T the elastic sig- 



nal decreases and eventually vanishes near Tc similar to 
an equivalent disorder signal in Zn-free LSC0.[l6|, |21| 
These results are consistent with a similar slowing down 
and subsequent freezing seen by /LtSR in underdoped 
cuprates j22l|. and agrees with the overall picture from 
NMR studies that Zn ions not only suppress dSC but also 
strongly enhance SDW correlations [23J. Recent neutron 
scattering ofFdetwinned YBa2Cu306.6 with 2% Zn found 
induced short-range magnetic order and a redistribution 
of spectral weight from the resonance peak to uniaxial 
IC spin fluctuations at lower energies 24 1 . 

Theoretically, the hourglass spin dispersion has been 
studied within various stripe models 25|, |26| , but a quan- 
titative description including both spin and charge de- 
grees of freedom, domain formation, glassiness, and/or 
slow fluctuations is still lacking. Some of these issues 
have been addressed recently 261432 1. but modeling the 
dynamics of the neutron response in the presence of dis- 
order remains limited at present 26|, 130 . |33| . 

Here, we focus on the disorder-induced magnetic phase 
known to exist in a correlated dSC and calculate the dy- 
namical spin susceptibility. For a homogeneous dSC our 
model reduces to a system of Bogoliubov quasiparticles 
whose magnetic response depends crucially on the pres- 
ence of a dSC gap|3a|. In the following, we investigate 
the role of spatially inhomogeneous local moments in- 
duced in this quasiparticle system, i.e. our formalism 
includes both the spin susceptibility from the dSC con- 
densate and the local moments. Thus it goes beyond a 
spin-only approach by including also the charge carriers 
important for these metallic systems. This is consistent 
with recent neutron experiments providing evidence for 
a two-component scenario where intertwined re gion s of 
itinerant carriers and local moments coexist |3J, |35[- A 
main goal of the calculation is to elucidate the spin fluctu- 
ation spectrum from an SDW droplet phase and coinpare 
it to the disordered stripe phase. 



Model. The Hamiltonian is given by 
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where c]^ creates an electron on site i with spin a, and 
tij = {t,t'} denote the two nearest neighbor (NN) hop- 
ping integrals, Vi is an impurity potential from a set of 
N point-like scatterers of strength V = IQOt, ^ is the 
chemical potential and Aij is the d-wave pairing poten- 
tial between sites i and j. The amplitude of Ay- is set 
by the dSC coupling constant g. Below we fix the pa- 
rameters t' = — 0.35f, g=:0.3t leading to pairing ampli- 
tudes A ^ O.lt, and adjust fi to give a hole doping 
X = 1 — n ~ 0.13. We solve Eq.dT]) self-consistently on un- 
restricted 24 X 24 lattices by diagonalizing the associated 
BogoHubov-de Gennes (BdG) equations at T = 0[37|. 

The model given by Eq.([T]) has been used extensively 
to study the competition between bulk dSC and SDW 
phases, the origin of field-induced magnetization, and 
moment formation around nonmagnetic impurities [38|. 
In the case of many impurities, Eq.([T]) was used to model 
static disorder- induced magnetic droplets |27l l30l - [33| . and 
to explain how these may increase in volume fraction 
when moving to lower doping levels and eventually form 
a SDW long-range ordered phase. More recently Eq.([T]) 
extended to the vortex state was used to obtain a semi- 
quantitative description of the T-dependence of the elas- 
tic neutron response in underdoped LSCO[33|. 

The transverse bare spin susceptibility Xo^ {rt , tj , uj) = 
—ij^ dte^^* ([tjf (i),tTj(0)]), can be expressed in terms 
of the BdG eigenvalues En and eigenvectors Un,Vn as 
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Including the electronic interactions within RPA we 
find for the full susceptibility 

X^^(rl , f„ c) ^^ [1 - Uxl'-iu:)]',!,,, Xo''(ri,rj,^)- (5) 

Fourier transforming with respect to the relative coordi- 
nate r = ri— fj defines the spatially resolved momentum- 
dependent susceptibility x('f, -R, w) = ^-e"^''x(^, r, w). 
Averaging over the center of mass coordinate i? = (r^ + 
rj)/2, this expression gives the susceptibility x(g,a;) rel- 
evant for comparison with neutron measurements. 



Results. The Hamiltonian ([T|) supports both a 
correlation- and disorder-induced SDW phase. Specifi- 
cally, in the clean case above a critical repulsion Uc2 — 
2.75i a global stripe phase is the favorable state. Below 
Uc2 the ground state is a homogeneous dSC but nonmag- 
netic disorder may locally induce moments if C/ > Ud 
where Ud is another critical interaction strength|30l.l37|. 



In Fig. [T]we show the magnetization and spin suscepti- 
bility in the clean system for U < Uc2 [Fig. [Ua,b)] and 
U > Uc2 [Fig. [llc,d)]. Figure [IJa,b) is a homogeneous 
dSC whereas the stripe phase in Fig. [TJCjd) coexists with 
a stripe-modulated dSC order parameter. The spin sus- 
ceptibility in Fig. [TJa) is dominated by a spin-1 in-gap 
resonance mode whose very existence depends crucially 
on the presence of a dSC gap[3g. By contrast, the spin- 
wave branches in the ordered case in Fig. [Ijc) are inde- 
pendent of the presence of a coexisting dSC condensate. 

We now demonstrate how disorder masks the clear dis- 
tinction between the susceptibilities in Fig. [Ua,c). In 
fact by increasing the concentration of nonmagnetic scat- 
terers one may tune the spin response from that in Fig. 
[IJa) to[IJc). This is evident from Fig. [5] showing the 
evolution of the spectrum in Fig. [ija) as a function of 
increased disorder concentrations rii. One clearly sees 
the impurity-induced low-energy spectral weight enter- 
ing the legs of the lower part of the hourglass, as well 
as an overall shifting down and broadening of the disper- 
sion branch. The shift is seen more clearly in Fig. [3l^a-c), 
which show w-cuts at qic = (|, l)''' for various values of 
U and Hi. At lower U, impurity concentrations of just a 
few percent may strongly enhance the low-energy weight 
and generate an in-gap mode [Fig. |31^a,b)], in agreement 
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FIG. 1: Intensity plot of xil,^) along a cut q = (g, tt) (a,c) 
and magnetization (b,d) in the clean system with (without) 
stripes shown in (a,b) [U = 2.6f] ((c,d) [U = 2.9t]). The 
susceptibility in (c) is the sum of xi'l,''^!^) and x('''"i9i'^) 
representing domain averaging. 



with neutron scatterin g on Zn-substitutcd samples, and 
in applied fields [H, \m^ H IH]. For sufficiently large 
disorder concentrations rii and repulsion U, the static 
component dominates x(9/CjW), as shown in Fig. ^c). 

The results shown in Figs. [2E] arise from the impurity- 
induced magnetism as can be directly verified by the spa- 
tial information contained in ^(q', i?,a;). Figures [3l^e,f) 
compare the dynamical susceptibility of the clean system 
x(g/c,a;) [Fig. ^c)] to the same quantity measured at 
the nearest neighbor site of a single nonmagnetic impu- 
rity [Fig. ^i)]. The local antifcrromagnctic (AF) droplet 
induced by the impurity is shown in Fig. |3|d). From Fig. 
[3ljf ) it is evident that the NN-sites to the impurity exhibit 
a low-energy peak which is responsible for the low-energy 
weight in the spatially averaged xiliCi'^)- 

Figures [IKS] apply to the so-called droplet phase where 
the system spontaneously generate impurity-induced lo- 
cal AF puddles as shown on the right panels of Fig. [21 
Note that the susceptibility presented here is for a few 
percent strong scattcrers representing e.g. Zn but that 
the overall phenomenology may also apply to intrinsically 
disordered cuprates like LSCO and BSCCOJlSlai. 

In the remainder we turn to the limit where U > Uc2 
and compare the results obtained above to a disordered 
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FIG. 2: (a,c,e) Same as in Fig. [TJa) with the addition of 
disorder. (b,d,f) show the associated real-space static mag- 
netization. Parameters are: U — 2.6t and rii ~ 1% (a,b), 
Tii ~ 2% (c,d), and rii — 3% (e,f). The results in the left 
panels are averaged over ten random disorder configurations. 
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FIG. 3: (a-c) x{qic,^) as a function of uj for difi^erent rn 
and U = 2.5t (a), U = 2.6t (b), and U = 2.7i (c). (d) 
shows the real-space magnetization induced by a single non- 
magnetic impurity at R=(12,12) for U — 2.7t. (e,f) display 
the susceptibility at the IC wavevector qic = (f , 1)'"" in the 
bare {U = 0) case (blue solid curve), and the full local suscep- 
tibility at _R = (11, 12) (red dashed line) with U = 2.5t < Ud 
(e) and U — 2.7t > Ud (f). The black solid lines in (e,f) show 
the spatially averaged full RPA susceptibility x(<J/c,'^)- 



stripe phase |26l433l |. Figure |4l^a,b) show how the mag- 
netization from Fig. [TJd) evolves into a nematic phase 
when disorder destroys the translational invariance. For 
a discussion of the disorder effects on the spin spectrum it 
is convenient to look separately at the dispersion perpen- 
dicular [Fig. |4ljc,d)] and parallel [Fig. |4l[e,f)] to the stripe 
direction. In the limit of low disorder [Fig. |ll[a,c,e)] the 
spin fluctuations are qualitatively similar to the clean 
stripe case from Fig. [IJc) (which shows the directional- 
averaged susceptibility). When sufficient disorder is 
added, the stripes break up and meander strongly, which 
is revealed in the spin response as a significant broaden- 
ing of the high-energy part of the hourglass including the 
spin resonance at (7r,7r), and a shifting-down of spectral 
weight in the response parallel to the stripe direction as 
seen from Fig. IH^f). A natural question arising from the 
similarity between the spin spectrum in the droplet phase 
[see e.g. Fig. EJ^e)] and the disordered stripe phase in Fig. 
|4] is how to experimentally distinguish between them. 
Since the droplet phase preserves the C4 symmetry of the 
underlying lattice whereas the stripe phase breaks C4 to 
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FIG. 4: (a,b) Real-space magnetization witii 2% (a) and 4% 
(b) disorder added to a stripe ordered system similar to Fig. 
[TJd). (c-f) show the associated susceptibility transverse (c,d) 
and along (e,f) the stripe direction for the same impurity con- 
centrations, i.e. 2% disorder in (c,e) and 4% in (d,f). 



C2 which remains robust even with significant disorder, 
we find that q^ — Qy anisotropy in the neutron response 
is the best way to answer this question. Experiments on 
detwinned underdoped YBCO samples have recently re- 
ported such anisotropy in the low-energy intrinsic |3| and 
Zn- induced spectral weight j24| . Hopefully future exper- 
iments may reveal whether similar results apply to e.g 
detwinned LSCO crystals. Careful studies of any poten- 
tial anisotropy in the inelastic signal may also reveal to 
what extent the droplet pha se can be viewed as a pre- 
cursor of the stripe phase [33[. 

Conclusions. We have studied disorder effects on 
the hourglass spin fluctuation spectrum relevant for the 
cuprates. We have focussed on two regimes of electronic 
repulsion U leading to a magnetic droplet phase for low 
U and a disordered stripe phase at larger U. Within 
this scenario, it is likely that the underdoped state is 
well described by our results for larger U, since we ex- 
pect that the effective interaction scale U/t increases as 
one approaches the Mott limit. Wc have also shown how 
neutron studies of detwinned samples may be used to dis- 
tinguish between them despite their surprisingly similar 
fluctuation spectrum in as-grown samples. 
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